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Abstract: We consider a physically viable cosmological model that has a field dependent 
Gauss-Bonnet coupling in its effective action, in addition to a standard scalar field 
potential. The presence of such terms in the four dimensional effective action gives rise 
to several novel effects, such as a four dimensional flat Friedmann-Robertson-Walker 
universe undergoing a cosmic inflation at the early epoch, as well as a cosmic accel¬ 
eration at late times. The model predicts, during inflation, spectra of both density 
perturbations and gravitational waves that may fall well within the experimental 
bounds. Furthermore, this model provides a mechanism for reheating of the early 
universe, which is similar to a model with some friction terms added to the equation of 
motion of the scalar field, which can imitate energy transfer from the scalar field to matter. 
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1. Introduction 

Although Einstein’s theory has been proven to be remarkably simple and successful as a 
classical theory of gravitational interactions, there are several observational facts which it 
has failed to elucidate. These cosmological conundrums include both cosmic inflation, or 
a period of accelerated expansion in the early universe, and a recent acceleration in the 
expansion of the universe. 
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Inflation in the early universe is a very attractive proposal for explaining the present 
large scale homogeneity and high degree of isotropy of the universe (one part in 100,000), 
in addition to the observed spectrum of density perturbations, which is usually attributed 
to a scalar field rolling down a shallow potential. Similarly, the current acceleration of the 
universe, as indicated by recent cosmological results [1], is usually attributed to some form 
of cosmic fluid having a large and smoothly distributed negative pressure, usually called 
dark energy or dark pressure. 

Cosmologists have long wondered why/how the universe has been recently accelerating: 
is it due to a pure cosmological constant term, or due to some sort of negative pressure gen¬ 
erated by one or more dynamical scalar fields, or something else? In recent years, different 
explanations have been provided for both inflation and the current epoch of acceleration: 
some examples of recent interest include brane-world modification of Einstein’s general rel¬ 
ativity (GR), including a 5d DGP (Dvali-Gabadadze-Porrati) model [2]. The names of the 
dark energy candidates run the gamut from f{R) gravity [3] (modifying in a very radical 
manner the Einstein’s GR itself) to ghost condensates (the idea which ncludes a more or 
less disguised non-locality) [4]. Many of these new proposals are pathological and do not 
appear more appealing than the two long envisioned alternative models of dark energy: a 
cosmological constant [5] and a slowly varying A-term [6]. 

The cosmological constant is a pure dark energy or vacuum energy, while the variable 
A-term is some kind of exotic matter or a slowly varying potential of a scalar field, usually 
referred to as quintessenee [7]. This last category can comprise a Gasimir energy or vacuum 
polarization effect from additional compact or curved non-compact spatial dimensions, 
which only weakly couple to ordinary matter in contrast to most tentative quintessence 
models, including k-essence [8] or curvature quintessence. In a brane world scenario [9], for 
instance, the vacuum energy (or the dark energy) may be viewed as a smooth brane-tension 
if our universe is a 3-brane embedded in higher dimensional spacetimes. 

Indeed, the past decade has witnessed signihcant progress in the building of inflationary 
models as extensions of standard cosmology, to accommodate the effects of dark energy. 
However, most of the inflationary type potentials studied in the literature are picked up in 
a very ad hoc fashion, rather than constructing such a potential as a valid solution of the 
held equations that follows, for instance, from low energy string effective actions. In this 
paper we initiate work in this direction. 


2. Action and equations of motion 

For a given scalar held a with a self-interaction potential V (a), an effective action as the 
low energy approximation of a fundamental theory of gravity and helds is written as 


s= d^x 
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( 2 . 1 ) 



where k is the inverse Planck scale Mp^ = and cr is a classical scalar field 

whose stress-energy tensor acts like a time-varying A. When studying the dynamics of 
an inflationary universe, the choice of the field potential V{a) is of particular interest. 
However, the origin and exact nature of the field a and the functional form of V(a) that 
acts as an extra source of gravitational repulsion (or dark energy) are not precisely known. 
A real motivation for the gravitational action of the form ( |2.lD arises from the following 
interesting observation. Inflation with the dynamics of a scalar field with a self-interaction 
potential, V{a), provides a negative pressure to drive the accelerating expansion of the 
universe [10] as well as a mechanism for generating the observed density perturbations [11]. 

The effective action (^) is found to be remarkably simple, but it excludes, at least, 
one important piece, which is the coupling between the scalar field a and the Riemann 
curvature tensor. This is because the field cr, if its vacuum expectation value is to describe 
the size and shape of the internal compactification manifold, generically couples with the 
curvature squared terms [12,13]. To this end, one welcomes the idea that a dark energy and 
its associated cosmic acceleration is due to a modification of general relativity such that 
the scalar field a couples to gravity via the curvature squared terms in the Gauss-Bonnet 
(GB) combination. The effective action for the system may be taken to be 


5 = d^x 




( 2 . 2 ) 


where Q = R'^ — is the Gauss-Bonnet invariant and V{a) and ^(cr) 

are general functions of cr. Dots represent other possible contribution to the gravitational 
action, such as C(<7)(Vc7)'^ [13,14], which we drop here for simplicity. This is justified, since 
the quartic term (Vcr)^ usually decays faster than the GB curvature invariant. The most 
desirable feature of GB type curvature corrections is that only the terms which are the 
second derivatives of the metric (or their product) appear in the field equations: a feature 
perhaps most important in order to make a theory of scalar-tensor gravity ghost free. 

For some, a weakness of this approach to model dark energy may reside in the standard 
argumentation for this particular kind of modification of Einstein’s theory, especially when 
one remembers that string/M theory predicts not only the fourth derivative gravitational 
term, like a GB invariant, but also higher-order terms. As is known, the field dependence 
of the coupling ^(u) has its origin in the variation of the background spacetime, and, in 
a spatially flat spacetime background, the GB term is subject to a non-renormalization 
theorem which implies that all moduli dependent higher loop contributions (e.g. terms 
cubic and higher order in Riemann tensor) vanish [12,15]. 

We anticipate that the Gauss-Bonnet invariant Q decays faster than the coupling ^(cj) 
grows, so that the term ^{cr)0 is only subdominant in the effective action. Moreover, the 
Gauss-Bonnet term is a topological invariant in four dimensions if ^{a) is a constant, but 
not if ^(cr) is a dynamical variable. 
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In an influential paper, Antoniadis et al. [15] demonstrated the existence of cosmolog¬ 
ical solutions which avoid the initial singularity and are consistent with the perturbative 
treatments of the string effective actions; see, e.g.; [16] for other generalizations. The es¬ 
sential ingredient of their method is a field-dependent Gauss-Bonnet coupling, ^(cr), where 
a characterizes the overall size and the shape of the internal compactihcation manifold. In 
the present paper we go one step further by demonstrating with new exact non-singular 
solutions that inflation in the early universe, as well as a cosmic acceleration at late times, 
may be explained by introducing a scalar potential V (a) for the modulus field a, in addi¬ 
tion to a field dependent Gauss-Bonnet coupling, ^(cr). We have reported some of the key 
results of our work in ref. [17]. The present paper will further elaborate the details of the 
model. For mathematical simplicity, we henceforth define the function /(cr) = A — ^(cr). 

While the field potential, V{a), was absent in some original string amplitude compu¬ 
tations, e.g. [12,13], implying that V{a) is a phenomenologically motivated field potential, 
it is quite possible that, in the presence of additional sources (like branes, fluxes), the 
string effective action would incorporate a non-trivial field potential, F(cj), as is the case 
revealed recently from string theory cosmic landscape [18,19]. Unfortunately, we do not 
have a precise knowledge about the field potential V{a). In string theory context, any 
such potential might take into account some non-perturbative effects, as arising from the 
effects of branes/fluxes present in the extra dimensions. In this sense, one may consider 
our model as string-inspired. 

An interesting question that we would like to ask is what new features would a dy¬ 
namical Gauss-Bonnet coupling introduce, and how can it influence cosmological evolution? 
Recently, it was suggested in [20] that the action (|2.2D with some specific choices of the 
potential and the scalar-GB coupling, namely, 

U((T) = Uoe-"(*)/"°, /(a) = /oe"W/'^°, (2.3) 


may be used, in four dimensions, to explain the current acceleration of the universe, in¬ 
cluding the phantom crossing the phantom divide {w = —1), with effective (cosmological 
constant or quintessence) equation of state of our universe; also see [21] for other interesting 
generalizations. 

In this paper, instead of choosing particular functional forms for V{a) and /(cr), as 
above, we present exact cosmological solutions that respect the symmetry of the field 
equations which follow from (|2.2|). The graviton equation of motion derived from the 


action (^) may be expressed in the following form (see, for example, ref. [22]) 


+ \9^.uV{a) + 2{g^,V^ - V^V,)(/(ct)R) - 4g^,V^VP{f{a)Rxp) - 4V\f{a)Rp,) 

+ AVPVp{f{a)R,p) + AVPV,{f{a)Rpp)+AV^f’V^\f{a)Rpp,x), (2.4) 
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where = RR^y + Rp,paxRv — ^Rp. ^Rvp- The equation of motion for the scalar field 
a is similarly given by 


0 = yV^fj 


dy{<^) . df{a) 

da da 


(2.5) 


Next, we consider a four-dimensional background spacetime defined by the standard Friedmann- 
Robertson-Walker metric: 

3 

= —dt^ -|- a{t)‘^ ( 2 . 6 ) 

i=l 

In this background, 

g = 2AH^{H + H^), (2.7) 


where H = d/a is the Hubble parameter and d = ^. The quantity 4:Xpy — gpyG vanishes. 
The (/i, u) = (t, t) and (x, x) components of the field equations have the following forms 


Q + 'la'^+ V{a)-2AaH^^, (2.8) 

K A (L(7 

0 = ^ (2iJ + 3if=) + (4 + <^^4) + + H-) + - V(.). 

(2.9) 


The time evolution equation for a{t) (cf equation (plq)) can be written as 




dt 


dt 


7 (fjfj -b 3idcT^) -b ^ {V{a) - f{a)g) + = 0 

^ g + A(a)) + 6id + d = 0, 


where we have defined 


k{a) = V{a) - f{a)g, 6 = f{a) 


dt 


( 2 . 10 ) 


( 2 . 11 ) 


We will call A(cj) an effective potential. Due to the Bianchi identity, one of the field 
equations (|2.8|) - (|2.10|) is redundant and hence may be discarded without loss of generality. 
In the limit /(cr)id^ —> 0, the action (^) reduces to (2.1). 


The 5 term in (p.lC ) may account for the creation of particles due to time variation of 
g. A friction-like term like this was first introduced in [23] on phenomenological grounds 
and a question was subsequently raised in [24] about the physical origin of such term; in 
our model, the 5 term is a clear manifestation of a non-trivial coupling (or back-reaction) 
between the field a and the time-varying Gauss-Bonnet curvature invariant. This term 
represents a clear advantage of (|2.2D over (^) as a cosmological action. 
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3. Scalar field as a perfect fiuid 


In a spatially flat background, the Einstein tensor, has components Goo = 3JI^, Gu = 
—a?‘{2H + 3//^). Assuming that the stress-energy tensor T^i, is described by a perfect fluid 
of the form Too = p, Ta = we hnd 

2H + 


P = 


,2 ’ 


p = 


We also immediately hnd from equations 


•(H) 


p = 1^2 + E(a) - 

p=\d’^- V{cj) + 


(3.1) 

(3.2) 

(3.3) 


where / = /(u). Note that this differs from the usual expression, due to the held-dependent 
coupling constant, /(u). To guarantee that the energy density of the scalar held is always 
positive, we require fH < 0 . 

We assume the scalar held obeys an equation of state, with an equation of state pa¬ 
rameter given by 

p 2, 2g-l , , 

w = - = -l--h = —-—, (3.4) 

p 3 3 

where q = —aaja^ is the deceleration parameter. A minimal criterion for getting an 
accelerating expansion is p -|- 3p < 0 , which is clearly a violation of the strong energy 
condition for some time-like vectors the latter states that or 

equivalently p -|- 3p > 0 and p -|- p > 0. Note that 


p + 3p = 27d^ - 2V{a) + 2A-{NH^f), 


(3.5) 


where N = f H dt, so that dN/dt = H. For / = 0, the acceleration condition p -|- 3p < 0 
holds when V{a) > 717 ^. In the case / 7 ^ 0, however, whether the condition V{a) > 
is sufficient or not depends on the sign of the time-derivative of the coupling, /(u). For a 
canonical scalar held (i.e. 7 > 0), both / < 0 and H < 0 hold, in general. The acceleration 
condition p -|- 3p < 0 is satished for f{a) < 0. Similarly, we hnd, 

d 


P + P = 73-^-|-8 




(3.6) 


and not zero as it would be in a spacetime which is exactly de Sitter: 7 = 0, /(a) = const. 
The null energy condition ^ 0 (for some null vectors x^)j or equivalently p + p > 0 

may be violated by allowing 7 < 0, or instead by taking ) + f jd^dt < N. 


4. General solutions 

To simplify the study of the model, we dehne the following dimensionless variables: 


X = 


7«: 


y = 


K^V (a) 


u = , h = 


H 

IP' 


(4.1) 
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The equations of motion, (^■^)- (|2.10 ), then form a set of second order differential equations 
(see the Appendix). In the absence of the GB coupling, so u = 0, we find a simple 
relationship; 

y = 2) + h, X = —h. (4-2) 


A physically more intriguing case is u 7 ^ 0, for which the variables satisfy 

y = 3 + h + - \u' + (5 — h)u — 2{h^ + h' + 5/i)tt] , (4-3) 

X = —h — - \u' — (1 + h)u — 2{h? + h' — /i)u] , (4-4) 

where the prime denotes a derivative w.r.t logarithmic time or number of e-folds N = 
ln(o(t)/oo). Thus there can exist a large class of solutions with different u{N). 

The quantity A(ct) = V{a) — f{cr)Q may act as an effective potential for the field a. It 
is clear that A((t) is a second order differential equation with non-constant coefficients: 


K^A{a) = {y — 3m(1 -|- h)) 

= \u'' + — h)u — 2 ( 8/1 + + h' + 3)u] -|- + h). (4.5) 

When solving (^ ) for n, one finds that the homogeneous part of the solution corresponds 
to solving the complementary differential equation, 


u" -|- (5 — h)u — 2 ( 8/1 + h'^ + h' + 3)u = 0. 


(4.6) 


The homogeneous solution is the full solution when one makes the ansatz 

K^A{a) =H‘^{3 + h), (4.7) 

which removes all nonhomogeneous terms. Of course, /(u) = 0 (i.e. tt = 0) is the triv¬ 
ial solution for ( [4.6D , which corresponds to the absence of Gauss-Bonnet coupling. Note 
that by setting k^A = {3 + h) one effectively makes the ansatz, V{a) = f{cr)G -|- 

MpH‘^{a) {3 h), reducing the number of arbitrary parameters of the model to one; fixing 

h alone will fix the function u{a{N)), or vice versa. A salient feature of our construction of 
cosmological solutions is that, while the contributions coming from both the field potential 
V{a), and the GB potential term, Vgb{^), may be large separately, the effective potential, 
A(it), can be exponentially close to zero at late times, as it dynamically relaxes to a small 
value after a sufficiently large number of e-folds of expansion. 

A commonly discussed alternative is the following. One solves the (modified) Einstein 
equations by making assumptions about the (functional) form of the field potential V{a) as 
well as the GB coupling /(u), which may be motivated by the leading order terms obtained 
(see, for example [20,25]), by a time-dependent (cosmological) compactification of classical 
super gravities. However, in this fashion, one may find a good approximation at each energy 
scale, but the corresponding solution will have little relevance when one attempts to study 
a wider range of scales, including an inflationary era. 
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4.1 Homogenous solution with constant h 

As a first reasonable approximation at low energy, we make the ansatz h = H/H^ = 
H'/H ~ const < 0. More specifically, when h = ho, the Hubble parameter is given by 

H = = {co-hot)-\ (4.8) 


where Hq (or cq) is an integration constant. Clearly, for t > 0, we require ho < 0 (and also 
Hq > 0). The scale factor is given by a{t) = ao(co — hot)~^/^°, implying that the universe 
accelerates when — 1 < /iq < 0. For t ~ 0 (and/or /iq — 0), the size of the universe in this 
regime grows approximately as e^*, where H ~ Hq ~ Cg Imposing (|4.6| ), we find the 
homogenous solution: 

u = ui e"i^ + U 2 (4.9) 

y = 3 + /iQ + 3(1 + ho') {ui + U2 , (4.10) 

X = — — ^7 + 5ho + \J 9/ig + 54/io + 49^ ui 

- ^ (7 + 5/io - y^9/i2 + 54 ;,^ + 49 ^ - ho , (4.11) 

where 

Cl = — — 5 + 9 /iq + 54/10 + 49^ , 

C2 = — Tho — 5 — ^9 /iq + 54/10 + 49 J , (4-12) 


and ui,U 2 are arbitrary at this stage. Whether the coupling function u is decreasing or 
increasing with N depends on the value of ho- In particular, for ho ~ 0, we have ci ~ —6 
and C 2 ~ 1, in which case the second term on the rhs of ( |4.9| ) increases with N. To 
rescue it, one might require to set U 2 ~ 0. However, for /iq ~ — 1 (and hence w r-^ — 1/3), 
since ci —> —2 and C 2 ^ —4, only the first term on rhs. of (|4.9f ) is relevant for large N, 
unless U 2 ^ ui- For purely exponential solutions we require that ui and U 2 be real and 
9/10 + 54/10 + 49 > 0. This means that either 


ho < 



-4.88 


or 


ho ^ —3 + 


4v/2 

~ 3 ~ 


- 1 . 11 , 


(4.13) 


implying that w > 2.2 or tc < —0.26. However, only for ho > —1 (and hence w < —1/3), 
the universe is accelerating. 

On the other hand, if 9 /iq + 54/io + 49 < 0, then the solution is written as 


u{N) = (^iVY/-(9/i2 + 54/10 + 49) 

+ k 2 cos I N■\J—{9hQ + 54/10 + 49) 


(4.14) 
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where ki, k 2 are real. We see that in the latter case for /iq < 5 the GB contribution at late 
times (large N) becomes vanishingly small. This is automatically satisfied as 


-4.88 < /in < -1.11, 


(4.15) 


implying that —0.26 < w < 2.2. Equations (^ ) -( O^ provide, in some sense, a natural 
solution for Einstein’s equations, as the effective cosmological constant, A(cr), is dynamical 
and becomes arbitrarily small with time, while simultaneously they provide the simplest 
forms for the GB scalar coupling, /(cr), the held potential, V{a), and the kinetic energy, 
K(a). However, this scenario cannot hope to effectively describe more than one epoch in 
the history of the universe, as h is not dynamical. 


4.2 Homogenous solution with dynamical h 

Motivated by the form of the previous solution we now make no assumptions about the 
form of the (slow-roll type) variable h, but instead consider the ansatz that during a given 
epoch we may make the approximation: 

u{a) = f{a) H^ = {X- C{a))H^ ^ ut , (4.16) 

where the time, t, can be early or late, where \ueariy\ \'U'iate\ and loeaWyl > Wiate\- 
The Gauss-Bonnet invariant, Q = + H) = 24 which is positive for the 

accelerating solutions (i.e. o > 0), decays faster than the coupling ^{a) blows up, but it 
may well be that decays slower than /(u) (or ■^(cr)). In the case /(cr) , a ~ 0. 

One may motivate the above ansatz for the coupling u{a) in two different ways. Eirstly, 
the function /(cr) as a solution to a second order differential equation in u{N) can have two 
different branches, which may dominate at different time scales, as can be seen above in 
the case h = H ~ const. Secondly, in a known form for the coupling /(cr) (= A —^(cr)), 
derived from heterotic string theory, the function ^(u) may be given by [15,16] 

,^(cr) = (51n[2e‘^r/‘^(ie°^)], (4.17) 

where h is a constant and i/(6), with b = is the Dedekind //-function which is defined 
by t/(6 ) = (l “ As can be seen from the plots in figure |I[ ^((t)/( 5 can 

be well approximated by ln2 — ^cosh((T), while {d^/(ia)/5 by —^sinh|cr|. 

In the expansion of the function ^(cr) (and hence /(cr)), the term proportional to 
can dominate at early times, cr <C 0, while the term proportional to e°^ can dominate at 
late times, u S> 0. The fact that such an expansion is typical of string effective actions 
represents, in our opinion, an interesting aspect of such models and implies the existence 
of two periods of accelerating expansion of the universe. Though the function ^(cr) is 
symmetric about a —> —a, the coupling u{a) is not, since the Hubble parameter H is a 
monotonically decreasing function of N (= ln(a)), or the field cr, and H^ariy lO^^eV > 
~ 10-33eV. 
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da ' 





0 


Figure 1: The function ^(u) (solid line) is symmetric about a —a, while its first derivative, 
d^/da (dotted line) is antisymmetric. Note that in these plots we have set J = 1. 


By solving equation (4.7), or equivalently equation (4.6), we find 


h{N) = —(3 + /3tanh/3AA^, 


(4.18) 


where /S.N = N — Nt, with Nt being an integration constant, analogous to In(ao) in (A.5), 
and 

16 + a ^ V9a2 + 72a + 208 


/3 = 


P = 


(4.19) 


4 ’ ' 4 

Here we have suppressed the subscripts referring to the time on a. Because of our ansatz 
for u{a{N)) of the form (4.16) the solution (4.181) is a good approximation at early or late 


epochs only. Of course, one could make a more complicated ansatz for the function u{a), 
motivated by some specific particle physics models, and find the corresponding solution for 
h{N). We believe that the solution we have found above is sufficiently simple to explain 
the evolution of the universe both at early time N > N^^ariy ^-nd at late times N > Niate- 


From (4.18), we easily find that 


H{N) = ef = Ho e“ cosh (SAN, 


(4.20) 


where, again, the time t can be early or late. Writing the expressions for y{N) and x{N) is 
straightforward, using eqs. (Ei)-ii- Analogous to an inflationary type solution induced 
by a conformal-anomaly [26], the solutions given above are singularity-free. The effective 
potential therefore takes the form 


H 


2 r 




A{a{N)) = -^ 3-/3 + /3tanh/3AiV {coshpANf . 


(4.21) 


This possesses a local maximum around 0.15 < AN < 1, depending upon the value of the 
slope parameter, a (see figure ^ and generically decreases to zero with increasing N. In 
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fact, A(c7 ) is almost flat for AN > 1. As we are assuming that Nearly ^ A'late, the local 
maxima and ‘global minima’ of the effective potential can be at different energy scales 
for the early-time and late-time universes. The kinetic energy, K{a), decreases with AN 
for 6 < a < 0, and increases for 0 < a < 1. The dimensionless variable, x = K/H"^, is 
approximately constant with AN > 10 for a wide range of a. 



-6 "5 -4 -3 - 2-10 1 


Figure 2: The contour plots of the effective potential A{(7)/Hq with the height in log scale. The 
single solid (blue) line denotes where dA/dN = 0, giving the local maximum of the potential w.r.t. 
AN. For larger AN, the potential generically decreases exponentially to zero. 



Figure 3: The kinetic term K/Hq (left plot) and the dimensionless variable x (right plot), in 
logarithmic scales, as functions of a and AN, where x = {'^/2){&/H)^ = K/H^ = (7/2 )(t'^. 
For a > 1 (or a < —6), K < 0 and hence a; < 0, leading to a phantom type cosmology. For 
—6 < Of < 0.45, K/Hq rapidly approaches zero. In all plots, where applicable, we have set Ut = 10. 


One may express the field potential V{cr) and the coupling constant /(cr) as functions 
















of the field a itself. As the second plot in figure ^ shows, for a < 0.45, the parameter 
x{N) is almost independent of N, implying that N = Ao(T+ const; we will choose this last 
constant to be A*. For AN ^ 0, Aq (= ^ function of the slope parameter 

a only. This leads to the following expressions for the potentials; 


V(a) = Vbe“2^^°'"(cosh(/3AocT))^ 


f(<^) = 


X 

^0 


3 — /3 + /? tanh(/3Ao(T) + 3ut + /3 tanh(/3Ao(T 

ut e(2/3+“)Ao^ (sech(/3AofT))^ e^^^*, 


,^Q:(Ao(t) 


(4.22) 

(4.23) 


where Vq = MpH^ e and ut = ut e"'^*. 



Figure 4: The ratio A{a)/2K{a) in logarithmic scale, in different ranges for a. 


For —6 < a < 1, so /3 > j3, the Hubble parameter is a smoothly decreasing function of 
N in this regime. Typically, for a ~ 1 (or a ~ —6), since /3 ~ /3, we have H ~ Hq, /i ~ 0 and 
hence w ~ —1: a stage where A(cr) acts as a cosmological constant. While, for 1 > a > —2 
(or —6 < a < —4), we have 0 > h > —1 and hence —1 < w < —1/3. When the value of a is 
decreased from unity, the kinetic term lowers towards zero at large N, while the parameter 
x{N) is almost constant (cf figure For a < 0.28, the condition A(it) > 2K(a) holds in 
general. But in our model this does not necessarily mean the existence of an accelerating 
phase; as the figure ^ shows, A(fT) > 2K{a) holds even if —2 > a > —4; in this last case 
we get h < —1 and hence w > —1/3, implying a non-accelerating universe. 

We will consider two epochs, an initial inflationary epoch where N is assumed to grow 
from an initial value Aj < Nearly, and a late-time deceleration/acceleration phase where 
A becomes comparable to Niate, and A > Niate eventually holds. The universe starts to 
inflate when A > Nearly + 0-5. Consequently, at A > Nearly, we have a stage of inflation. 
For P{N — Nearly) ^ 2.5, the scalar field begins to freeze in, such that w < —1/3; the actual 
value of w depends on the value of slope parameter a, see [17]. After a certain number of 
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e-folds, AN = Nf — Ni, our approximation, ( 4.161) , with the time being early breaks down. 
For N < Niate, the universe is in a deceleration phase which implies that inflation must 
have stopped during the intermediate epoch. Sometime later, when N > Niate, subsequent 
evolution will be controlled by (|4.16| ) with the time being late. 

For PAN > 2.5, the variation in h{N) is small, so we can approximate it by 


h 


-P + P (A1V>0) 
-P-P (A1V<C0). 


(4.24) 


Rearranging this approximation for h, we find a = {ci, C 2 }, as defined in (4.12). This tells 
us that the earlier approximation, /i ~ /iq, is a good ‘early time’ or ‘late time’ approximation 
during a given epoch. Also note that for h < 0 (and hence w > —1) we require —6 < a < 1. 
Picking a. outside this range leads to ‘big rip’ type cosmologies for which h > 0 (and hence 
w < —1) in some regions of field space. 

4.3 Relaxation of dark energy 

There might be a large shift in the Hubble expansion rate during inflation, viz Hfee/ore ~ 
10^3 eP > Hafter 10 33 eP. One could therefore ask whether the solutions we presented 
above explain a dynamical relaxation of vacuum energy (or scalar potential) to the present 
value of dark energy, 10 Mp, after a significant period of inflation. This is quite 
plausible in our model. To quantify this, one considers the ratio of the Hubble parameters 
before and after N e-folds of inflation, which is given by 


e = 


cosh/3(AA-|-A) e 
cosh/3(AA) 


(4.25) 


As an illustration consider that a ~ 0.0143 ~ 1/70, so that qini — 0 at AN = Nini — No ~ 
0.33. Assuming 70 e-folds of inflation, i.e. N = 70, we find 


1.36 X 10 


-12 


(4.26) 


This value represents a shift in Hubble expansion rate, only during an accelerating epoch 
for which the initial value of the deceleration parameter is zero. Indeed, something like 
70 e-folds of expansion, as usually considered, is the minimum for inflation, based on the 
assumption of a constant Hubble rate. Practically, one requires much more expansion 
than e™ between the Planck time and the present. As emphasized in [27], one might need 
N = In > 70, rather than N = In(ajaj) > 70, to solve the various cosmological 

conundrums, including the flatness and horizons problems. The difference N — N, which 
is non-zero and positive, as long as —6 < a < 1, would characterize the extra amount of 
expansion (e-folds) required by the decrease of H during inflation. 

In our model, the total number of e-folds N required to get a small ratio, like e ~ 
10 “^®, depends on the value of the expansion parameter a, which may be related to the 


- 13 - 








slope of the potential coming purely from the Gauss-Bonnet coupling, VcsCo") = /(o') Q = 
3u{a)H‘^{l + h). This value would be minimum, N ~ 125, for a > —2 (or a < —4), while 
it would be large, /S.N > 0(300), for a < 1 (or a > —6). In the former case, inflation could 
occur slowly since ^(^) ^ 0, while, in the later case, it might occur more rapidly since 
^(^) <C 0. For example, for a = 0, the effective potential A{a) ~ 10“®Mp decreases to 
the present value of dark energy, namely, Aq 10 Mp, when Ntotai ~ 326. 

Assuming the universe has undergone a sufficient number of e-folds of expansion, like 
A^totai ^ 125, we find an effective scalar potential that can dynamically relax its value to 
the observed value of the cosmological constant, such that the field a evolves towards its 
minimum and the equation of state parameter falls in the range —1 < w < —1/3. 


4.4 Late-time acceleration 


The deceleration parameter q, which may be parameterized as a function of red-shift factor 
z, is given by 


= -(1 + ^)- 

The value of the function, h, for the current epoch must be determined by observation. For 
example, if the observed value of the deceleration parameter q{z) corresponds to ~ —0.6, 
then 


Qobs — 1 ^obs 

- -1 + Plate - Plate - “0.6 , (4.27) 

where we have made the late time approximation for hobs- Solving hobs = —0.4 we find 

aiate - -0.0148 or aiate - -5.3851. (4.28) 

Out of the two possible values of the slope parameter a, giving rise to the same value of 
deceleration parameter q{z), we should prefer to take the smaller value of a, in terms of 
its absolute magnitude, so that the coupling u{a{N)) ~ f{a)H^ varies only slowly with 
the expansion of the universe. In fact, the function u{a) is a measure of the strength of 
scalar coupling with the curvature terms: u{a) is likely to depend upon the gauge coupling 
strength. And, in the presence of matter sources, like standard model particles, the running 
of the gauge coupling could be small if |q:| < 0{^). 

Given the definition of redshift 


1. Z — Ojnowj (4.29) 

we can rewrite it as 1 -|- z = allowing us to plot recent acceleration versus z. At late 

times, A^o = Aliate, we may observe a significant variation in h over the range of red-shift 
z < 1.5 (cf figure |^, in excellent agreement with observation [1]. 
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Figure 5: The solution is modelled such that ln(l + z) = AN = In One can see a change from 
deceleration to acceleration as the red-shift factor z decreases in the range 0.5 < 2 < 1.5. 

4.5 Scalar potential: leading order behavionr 

It might well be that, asymptotically (i.e. at late times), both the functions h{N) and 
u{a{N)) take (nearly) constant values. In this rather special limit 



(4.30) 


we find the following solution 



(4.31) 


(4.32) 


(4.33) 


where aoK = sj (mq ~ (1 + uo)/ho) /27. It is not difficult to see that, to leading order in a, 
f{a) (X e'^/°'°. That is, the ansdtze like V{a) = and f{a) = characterize 

only the late time attractors, for which both u{N) and h{N) behave merely as constants. 
Consider a more general case for which h{N) = H/H^ = fig < 0, but u{a{N)) = 


uq e“^, a / 0. In this case, a is related to N via 



(4.34) 


where 


Cl = 2rio(a — 2/io)(l — /ig — a)- 


(4.35) 
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Both signs of a may be allowed as long as ci > 4ho e holds. Indeed, h(N) = H'/H ~ 0 
is special case, for which 

=vf 

and hence A((t) = 3MpHQ = Aq. Note that, since the equation of state parameter w = — 1 , 
A(it) acts as a cosmological constant term in the regime /i ~ 0. In this case a can be real 
even if 7 < 0, provided that both a < 0 and tto > 0 hold simultaneously. Nevertheless, one 
should be more interested in a canonical scalar (7 > 0 ), in which case one now requires 
0 < a < I and uq > 0. As we will see shortly, the slope parameter a in this range possesses 
an interesting feature that the spectrum of scalar (density) perturbation during inflation 
is almost flat, giving a scale invariant power spectrum, Ug — 1. 


5. Inflation and cosmological perturbations 

It is generally believed that during inflation the inflaton and graviton field undergo quantum- 
mechanical fluctuations, leading to scalar (density) and tensor (gravity waves) fluctuations, 
which in turn would give rise to significant effects on the large-scale structure of the universe 
at the present epoch. The spectra of perturbations may provide a potentially powerful test 
of the inflationary hypothesis. One could ask whether the model presented here will add 
to the search for a conceivable physical basis for inflation. We believe this is possible. 


5.1 Slow roll variables 


In a standard scenario, one defines the slow roll variables, in terms of the field potential 

F(u): 


1 (V'\^ _ 1 fV" V'a"\ 


(5.1) 


As above, the prime denotes derivative with respect to N, not the field a. These definitions 
for slow roll parameters may be justified since during the inflationary epoch V(a) 3> Vgb(o') 
holds. Note that Vgb(cj) = f{a')Q = 3 uH‘^{ 1 + h) ^ 0 as h ^ —1. Because of the flatness 
of V{a), a grows very slowly and essentially all the inflation occurs when V(a} » Vg'b(cj). 
Slow-roll requires that |e„ | <C 1, |? 7 « | <C 1. 

To impose the conditions on slow roll variables in a physically motivated and model- 
independent way, following ref. [27], we may define them in terms of the Hubble parameter 
H{a) and its derivatives: 


en 


VH 


2 (HaV ^ 2 fH'y 


2 2 

fH" 

H' a" 

H ' 

[if 

H a' 


(5.2) 

(5.3) 
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where, as before, primes denote derivatives w.r.t. N. One also defines the following 
parameter, which is second order in slow-roll expansion: 




1 


2k? 




1/2 


= ( 


f) 


1/2 


(5.4) 


(The defined above is not to be confused with the coupling function ^{a) we defined 
before; here we are adopting the notations which are standard in literature and we will 
not refer to ^{a) in this section). In fact, these definitions of slow roll variables may be of 
wider applicability than those defined in terms of V{a), as they are based on the fact that 
inflation occurs as long as ;^(^) < 0 holds. That is, during inflation, the comoving Hubble 
radius, l/{aH), must necessarily decrease, so that physical scales can grow more rapidly 
than the Hubble radius. Inflation ends at ^ (|) = 0, where k = aH = aeH{a)e~^^ is, 
by definition, a scale matching condition, where Oe is the value of the scale factor at the 
end of inflation, at which a coming mode crossed outside the (cosmological) horizon. 

Let us analyze the last case considered above in some detail. First, note that, defining 
a —> we can always absorb the coupling constant 7 into the slow-roll parameters. 

Without loss of generality, henceforth we define e = eull-, 1 ^ Ih/i-, ^ ^ Ch/'I', similar 
arguments would apply to other quantities, like, scalar and tensor indices, defined below. 

At the level of approximation we are considering in this section, namely, V{a) /!> 
the tilt parameter, up to first order terms (e and 77 ), may be given by 


TT-s ~ 1 — 4e + 2rj. 


(5.5) 


(We will discuss below about the validity of such a relation in the presence of the Gauss- 
Bonnet term.) As is generally the case, a value < 1 is easier to produce than > 1 for 
the model we are considering here; ~ 1 is the value that makes the (scalar) perturbation 
small in all scales, see, e.g. [28], for a review. In our model, > 1 is possible only if one 
allows a > 1 (or < —6). This last demand is however not a physically motivated one since 
the kinetic energy of the field a is negative in at least some regions of field space. 

In our model, the slow roll variables e and r] (and also ^) can vary with both the 
exponent and coefficient of the coupling constant u{a{N)), i.e. a and ut, in addition to 
the number of e-folds, Nj — Ni = AN (cf figure ^). In turn, as the figures [^ and show, 
it is possible to get a value of close to unity in a wide range of ut (or AN) by suitably 
choosing a. Interestingly, as shown in figure |9|, a value of in the range [0.89,1] may 
easily be obtained by taking different combinations of a and ut- One obtains the value 
Us ~ 0.95, in excellent agreement with recent observational data from WMAP, by taking 


AN ~ 70, aearly ~ -0.01, Uearly ~ 22. (5.6) 

For a < 0, with a smaller value of AN, one also requires a smaller value of ut; when 
AN = 50 and a = —0.01, we find n* ~ 0.97 for u^ariy — 31, while n* ~ 0.95 for u^ariy — 18. 
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Figure 6 : The slow roll variables e and ry as functions of a and /S.N. For a small and positive 
a, the universe is decelerating before inflation, AA^ < 0 , but e (and also 77 ) quickly takes a small 
value (e, 77 ^ 1) for AA > 10 and a > —0.2, leading to inflation. In fact, e ^ 0 also for a ~ — 6 . 
Outside the range —6 < a < 1, we get e < 0, leading to a phantom type cosmology. The condition 
for acceleration to occur is e < 1, so the universe is not accelerating for —4 < a < —2. 



Figure 7: The spectral index Us as a function of slope parameter a and AN. Both for a > 0 and 
a <C 0, the spectral index is independent of AN. 


Of course, a positive value for the slope parameter, a, is also allowed, as long as a < 1 
holds. Note, for a > 0 and N > 0 the function u{a{N)) ~ e"^ grows with proper time, t, 
or logarithmic time N. But, in our construction, any contribution like this, coming from 
Gauss-Bonnet term is exactly cancelled with the homogeneous part (i.e. terms multiplied 
by u{a) and its derivatives) of the field potential V{a). Also, both the Hubble parameter 
H and u{a)H'^ can be slowly decreasing functions of the number of e-folds, N, or the field 
a. 

Note that if V(a) = 0, then Einstein gravity may not be an effective theory at low 
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Figure 8: The spectral index Hs as a function of a and ut, with a fixed value of AN = 70. The 
plots do not change much while including terms second order in slow-roll, viz, Tig — 1 — 4e + 277 — 
2(1 + c)e^ — ^(3 — 5c)e77 + ^(3 — c)^^, where c ~ 5.08, as defined in [27]. A small difference is that 
now a smaller value of Ut is required so as to get the same value of . 


U Ut 




Figure 9: The contour plots of Ut vs a in the range rig = [0.89,1], with a fixed value of AN = 70. 
The spectral index, rig, decreases (increases) from left to right in the left (right) plot; rig > 1 can 
be obtained only for a < —6 (or a > 1). For a > 0, more positive is the value of a, smaller will be 
the coefficient Ut giving rise to the same value of rig. For a << 0, however, rig is insensitive to the 
value of Ut- For a > 0, the signs of running of spectral index can be different between the Ut < 0 
and Ut > 0 cases. 


energy, since the term u{a)H‘^ can easily dominate the Einstein-Hilbert term, Rj k?, which 
is proportional to ‘iH^. Even though the time scale for such effect to occur can be extremely 
large, given that a is very close to zero (i.e. /(a) scales nearly as 1/H{(7)‘^) and the 

coefficient ut can be extremely small, this provides an additional justification for considering 
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an effective action with a non-trivial field potential, V{a). 


5.2 Generation of perturbations 

Note that for the solution (4.20) the scale factor of the universe evolves as 


(5.7) 


where m = (3 — (3 <1 and r is the conformal time. In the case of power-law inflation, such 
as this, the amplitudes for scalar and tensor fluctuations may be given by [29,30] 


A. = 


Ax = 


\h\ 


. X Ha' 

^(1 +0.46c-0.73,)^ 

2A^pH 

hy/iiMp 


(1 - 0.27e) 


(5.8) 

(5.9) 


Two remarks are in order. First, for each of these perturbations the value of H (and hence 
h) and a' (= ^) must be evaluated when the wavelength of the perturbation becomes of 
the order H~^. Practically, it is more convenient to specify them as functions of e-folds 
N before the end of inflation. Second, the above expressions, with Aq'^^ = aP = 1, best 
approximate the results in a model with a self interaction potential V(a) alone [29], i.e. 
without the GB coupling. In fact, in the f(a) ^ 0 case, ^o’s are generally functions of 
f{a) and 77, not unity. As long as the GB term is only sub-leading to the scalar potential 
V{a), the above results would be available to leading order. 

Following [32,33], one may calculate the coefficients related to scalar and tensor type 
perturbations 

.(T) _ 


where v = ^ + 


1—m 


and 

4i?2 


2h + 


\l Q{s) — 


2x + 


6F2 
1 + 2F 


F-f 
1 + 2F 
1/2 


2x + 


'1 + 2F 
1 + 3F 


6F2 
1 + 2F 


-1 


Q{T) — 1 + 2F, 

S(t)Q(t) = 1 + 2/, 


(5.11) 

(5.12) 

(5.13) 

(5.14) 


where F = K^fH and x = {.p)^■ l'^ limit T —> 0 (or / —> 0), one recovers the 

standard result for which S(s) = 'S(r) = 1) pOp = pi pQ(T) = 1- As can be seen 
at the level of the field equations, equations (|A.1|) - ([A.3| ), the condition 0 < jF] <C 1, or 
alternatively 

F = K^fH = — (^) = -{u — 2hu) <C 1, (5.15) 
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is equivalent to V{a) S> Vg_b(o'). This may be satisfied, for example, by choosing ut = 10, 
AN = 70 and a < 0. |T| < 1 implies that aH decays faster than the function ^ grows. 
Since / = —^ — d^, calculating the above quantities requires knowledge about 

how the scalar field a varies with time, which is model dependent. 

Since h < 0, a small negative value of F would normally decrease the value of 5 ( 5 ) 
but would increase the value of This is opposite for the tensor modes: 5 ( 7 ^) will 

increase but will decrease. Although Ag and both vary with AN, or the choice 

of scale k = aH at which a comoving mode crossed outside the horizon, the ratio At!A s 
is essentially independent of such a scale, or the number of e-folds, AN = Nf — Ni. The 
tensor spectral index may therefore be given by [30,31], 


i 2 r 


^ - 2 ^ 




(5.16) 


In the plots of figures 0 and one may find a significant deviation for F ~ —0.5. 
Interestingly, for the solutions that we found above, this last last condition is extremely 
difficult to satisfy since it requires u' — 2hu ~ —4, but u' — 2hu ~ 0_ at the end of inflation 
for a < 0 , and u' — 2hu > 0 for a > 0 . 

Let us discuss about the validity of the relation, Ug — 1 — —4e + 2r], in the presence of 
a non-trivial GB coupling. As shown in [34], in general terms, the spectral indices of the 
scalar and tensor-type perturbations may be given by 


Ug — 1 = 3 — 2ug, riT = 3 — 2vt. ( 5 - 17 ) 

In our case, since Vg = ^ + it is sufficiently clear that — 1 = —< 0. For the 
tree-level solutions found by Gasperini and Veneziano [35], the authors of [32] found the 
unpleasing result ~ 4 and ht — 3. The solutions given in [35] are characterized by the 
additional presence of a dilaton background with V{4>) = 0. This suggests, in our 
case, that the V{a) = 0 case is not illuminating at least in view of inflationary paradigm. 
Indeed, the solutions found in [35], namely, 

a(r) oc |rIi?!)(r) = —\/31n |r| (5.18) 


are non-accelerating, so it is inconsistent in using non-inflationary solutions to estimate the 
inflationary parameters, like Ug. 

In our model, it is quite possible that during an inflationary epoch the field a changes 
only slowly, so the quantity x{N) ~ = (u')^ remains almost constant over expo¬ 

nentially large range of wavelengths, leading to an almost flat spectrum of perturbations 
of metric. Fluctuation of the field a leads to the effect that the duration of inflation is 
increased due to a local delay of time near the exit from inflation, which may be given by 


5t 


da 1 
a 2'Ka' 


(5.19) 
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Figure 10: (a) (left plot) The logarithm of the scale Hq/Mp vs a, giving rise to the value As = 
2 X 10“®, which may be matched to the density contrast, Sh, at Hubble-radius crossing. A particular 
value of a fixes the energy scale Hq in terms of Planck mass, e.g., for a ~ 0.0284, Hq ~ 10“® Mp. 
(b) (right plot) The tensor-to-scalar ratio, ApjAs, as a function of a: from left to right Ut = 
30,20,10,5,1. This ratio is independent of the number of e-folds. The single (horizontal) line 
corresponds to Ut = 0, in which case Ap/As (and hence np) can be large (r < 0.43). 


AN 




Figure 11: The tensor index Ap as a function of a and m*, with fixed AiV = 70. The tensor modes 
are almost negligible, np ~ 0, for a > 0 (or a ~ —6), or equivalently, for e <C 1. 


This may lead to a local density increase such that [11] 

H 


Oh ^ — 
P 


27r a' 


(5.20) 


Inflation also leads to creation of pertnrbations of the field a with wavelength greater than 
. The average amplitnde of scalar perturbations during a typical time interval is 
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given by [36] 


(5.21) 


\5a{x) 

Indeed, this result, in its more rigorous form ^(cr^) = was independently obtained 
in [37]. A useful observational constraint is the following. If a changed very slowly during 
inflation, then H/a' remained almost constant over exponentially large range of wave¬ 
lengths. The cosmic microwave background constraint is such that [38] 

5h ~ 1.92 X 10“®. (5.22) 

If this quantity is to be matched (precisely) with the scalar-index Ag, as the results given 
in [28] suggest, then, in our model, we can reproduce 5h of this magnitude by taking, for 
instance, Hq r\j 10 and a ~ —0.01. 


6. Towards reheating in an inflationary universe 


It is essential to have a good reheating process after inflation in a cosmological model. 


Recall that in our model the time evolution of a is given by (cf equation ( 2.10 )) 


dt 


g.T= + AM) 


= -fifl g^^) - 


( 6 . 1 ) 


where A(cj) = V{a) — f{cr)G and <5 = /(c)^. Note that the effective potential for the 
dynamics of the scalar field is A(cr) rather than V(a). The first term on the r.h.s. represents 
the energy loss caused by the expansion of the universe and the 6 term represents the energy 
density per unit time which is drained from the field a though time-variation of the Gauss- 
Bonnet term. Physically, such a term is expected since the modulus field a, which obtains a 
non-zero mass due to its vacuum expectation value, (a), is coupled to the curvature tensor. 

Suppose that, initially, d ~ 0 and dA/dt ~ 0 (dA/dcr is not essentially zero there). 
The evolution equation o then implies that d ~ 0 and hence 




^4(No-N) 

1 + eKNo-N) ’ 


( 6 . 2 ) 


where Hq and Nq are some integration constants. Suppose we start inflation at a point 
where N < Nq ^ and h = H/H^ = H'/H ~ —1, so that the field equations satisfy 

3&H ~ — (, a ~ <C 2>Ha. (6-3) 

4 \ d(T / da^ 


For N > Nq, the universe experiences a nearly constant Hubble flow H ~ Hq, leading to an 
early inflation (exponential expansion) of the universe, i.e., a ~ The physical Hubble 

radius H~^ increases with the number of e-folds N as the field a rolls down its potential 

^One can define the scale factor as a = so that N = ln(a(t)) = uj{t) and No = ojo- 
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and gets out of its local extremum (a point of inflection or critical point). Inflation would 
be eternal (i.e. without an exit) if <5 ~ 0 holds for all times. However, this is generally not 
the case since the Hubble flow is damped by an adequate cosmic friction term, and most 
of the evolution of the universe would be described by the <5 7 ^ 0 or 7 ^ 0 solution. At 
this point, one also notes that the term <5 can change its sign between accelerating > 0 ) 
and decelerating < 0) solutions, since ^ = 24 That is, given that f{a) does not 

change its sign at the transition point, d = 0, the term 5 acts as a (positive) friction term 
during an accelerating phase > 0 ), while it is opposite during a decelerating phase. 

So far we have neglected the couplings of the scalar field to radiation (matter) fields. 
In the case the effective potential A(cj) possesses a local minimum, such couplings may 
cause the rapid oscillatory phase to produce particles, leading to reheating [23,24,39] 

To this end, we can introduce some matter helds, which constitute ordinary matter and 
radiation field released by the decay of the field a. For a homogeneous cosmological model 
with vanishing space curvature, we find 

(a) + 

4 (2^2 + am) =-6^ (2.^^)-(6.5) 

where pa = + V{(y) — The 5 term in ( |6.5| ) resembles a drag term, which 

transfers energy from the motion of a and dumps it in the form of a radiation background. 
The equation for the evolution of the energy density in radiation (particle) is given by 

Pr + 4:Hpr — 5 = 0 . ( 6 . 6 ) 

Once reheating is completed the universe enters a standard radiation dominated FRW 
phase: 

1 1 1 rn 

Po = oPO = nPr: Po ^ Hr -2 ~ T- V 

00 Qj a L 

The field a thereafter remains subdominant for most of the time and only at late times, 
N > Niate, when the potential becomes sufficiently shallow, does one get acceleration, 
w < —1/3. This would then make the model complete. 

A few remarks are in order. The first is related to a reheating process after inflation, 
the details of which depend upon a choice of the potential V{a) and the coupling /(cr). 
More generally, it is essential to know whether or not the effective potential A(cr) has a 
minimum. So far in this section we have not specified the form of A((t); our discussion of 
particle production above is rather qualitative. 

Note that for the solution given in subsection (4.2), A((t) does not have a local mini¬ 
mum, rather it has a local maximum (cf figure |T^). In this sense the conventional reheating 

^However, in non-oscillatory models, for example a model of quintessence characterized by a single 
exponential potential, the instant preheating mechanism proposed by Felder et al [40] might be more 
efficient for particle production, which works even for a potential without a minimum. 
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Figure 12: Approximate behaviours of the effective potential A(cr) (not to scale) at early (left plot) 
and late (right plot) times; we have chosen the free (integration) parameters -ffo.early and ilo.iate 
such that the plots mimic a realistic feature that A(cr)eariy ^ A(CT)iate- 


process may not work and an alternative reheating method need to be employed: an effi¬ 
cient method of reheating is the instant preheating proposed by Felder et al [40]. In this 
case the field cr decays when it rolls down the potential, thereby producing heavy particles. 
However, details of the reheating process may be somewhat different and complicated in 
our model for at least two reasons. Firstly, as compared to a standard model, there is an 
extra friction-like term, namely S, which is non-zero as long as 7 ^ 0. Secondly, in our 
ansatz for u(a) (= f{(T)H‘^), we assumed that u{(t) oc e“*^, where the time t can be 

early or late. In order to retrieve the full potential, we may have to sew these two poten¬ 
tials in some way, possibly creating a barrier (and hence a minimum) somewhere between 
them. Work in this direction is currently underway. Numerical studies show that for an 
ansatz of the type u{a) = Ueariy e"‘'^-|-uo(iV)-buiate where uq{N) ~ 0(1), A((t) can 

have an effective local minimum, where inflation could end. This is in accordance with our 
observation in subsection (4.2), that inflation must have stopped during the intermediate 
phase. In fact, some of our solutions presented in the Appendix ^ possess an effective local 
minimum. In any case, the instant preheating proposed in [40] is perhaps the most efficient 
method of particle production in our model, as it works even for a potential without a 
minimum. 

The second issue is related to the nucleosynthesis bound. As is known, for inflation 
driven by a scalar field, there exists a tight constraint on the allowed magnitude of = 
p < 0.1 —0.2 at the time of nucleosyntheis; see, e.g., [41]. This often places a constraint 
on the model parameters, for instance on the slope of the (effective) potential. One could 

®In the Appendix we do not make any specific ansatz for u{a) but instead allow some of the variables 
to take their limiting values, like x ~ const and h ~ const 
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ask whether or not the nucleosynthesis bound will not be violated in our model. To address 
this question properly, we would need to know the precise form of the (effective) potential. 

It is always possible to constrain some of the parameters in our model by allowing V(a) 
and /(cr) to take some specific (functional) form. In subsection (4.2), we found a solution by 
making a specific ansatz for /(cr) (or f{a)H‘^). For this solution, when the gravity is coupled 
to matter and radiation, the necleosynthesis bound, in terms of the allowed magnitude of 
Ho-, constrains the values of the parameters (3 and (3, namely 4|/3 — /3pAg > where 

n = 3 (4) for matter (radiation) and Aq = ~ W\f}^N> 2 - nucleosysthesis 

bound may not be violated for {(3 — /3)Ao > \/5, given that < 0.2; in arriving at 

this result we have assumed that the effective potential after a certain number of e-folds 
{f3AN > 2) is approximated by A((t) oc ^ since tanh(/3AA^) ^ 1. Also note that 

in this limit, V{a) oc with a different proportionality constant. 

An inflationary type potential that we arrived at from studying symmetries of the 
field equations is approximated by V{a) = Vq exp[— p(cj)k(t], with p being a function of 
the field cr or the number of e-folds N. Thus, the reheating process discussed in [42] may 
be useful in our model. Finally, we would like to note that our model might inherit some 
of the features of quintessential inflation proposed by Peebles and Vilenkin [43] in which 
the potential consists of two parts: X{a^ + M^) {a < 0) for inflation and \M^/{a^ + M^) 
{a > 0) for quintessence; both models can lead to tracker solutions at late times, though 
the forms of scalar potentials are quite different. 

7. Conclusions 

In this paper we presented an analysis of accelerating/inflationary cosmologies by introduc¬ 
ing in the effective action a field dependent Gauss-Bonnet coupling, other than a standard 
field potential for the field a. We find that the dark energy hypothesis fits into a low 
energy gravitational action where a scalar field is coupled to the curvature squared terms 
in Gauss-Bonnet combination. It is established that a GB scalar-coupling can play an 
important and interesting role in explaining both the early and late-time evolutions of the 
universe as well as providing a mechanism for reheating. 

That we are able to explain accelerating universes using exact cosmological solutions in 
a modified Gauss-Bonnet theory, leading to a small deviation from the w = —1 prediction 
of non-evolving dark energy (or a cosmological constant) is likely to have a serious impact 
in search of a viable dark energy model. Our work also provides extension of quintessence 
(or time-varying A) model in which part of the dark energy comes from a field dependent 
Gauss-Bonnet interaction term. 

One of the key results is this; in the absence of a Gauss-Bonnet coupling, the ten¬ 
sor/scalar ratio is usually non-zero. However, with a non-trivial scalar Gauss-Bonnet cou¬ 
pling, i.e., f{a) /, or effectively, u{a) = f{a)H^ ~ e“'^ 0, such a ratio can be negligibly 
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small if the expansion parameter a takes a small positive value, a > 0.1, and hence < 1 
and riT — 0, leading to Harrison-Zel’dovich spectrum. Unlike a naive expectation, the 
inclusion of a (scalar) field dependent Gauss-Bonnet coupling /(cr), in addition to a field 
potential V{a), into the effective action, could make the observability of tensor/scalar ratio 
and related inflationary parameters more achievable. 

We emphasize that, in contrast to previous analysis, our calculations have all been 
implemented by the functional forms of the scalar potential and GB scalar coupling, as 
suggested by the symmetry of the held equations, rather than choosing particular model 
dependant forms for them. We have given in the Appendix the exact solutions for some 
special cases, about which a general comparison can be made in terms of the homogeneous 
solutions we presented in the bulk part of the paper. 

Regardless of whether the model studied here appears natural or otherwise, it should 
be observation that determines whether or not it is correct. The current and future obser¬ 
vations might make stronger demand on theoretical precision of inhationary parameters, 
including, the scalar and tensor spectral indices, and are certain to constrain a number of 
parameters of our model tightly, including the Gauss-Bonnet coupling constants. 
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8. Appendix 

In terms of the dimensionless variables dehned in equation ([l.lj), the equations of motion, 
(2.5)-(2.1C), form a set of second order differential equations: 


0 = —3-|-x-|-y — 3{u — 2hu), 


where 


so that 


dX dX 
A = = 0,—r- 

dN da 


N = / Hdt = In 



(A.l) 

2h‘^u + X — y + 2h + 3, 

(A.2) 

h + l){u — 2hu). 

(A.3) 

1 dX 


hIF' 

(A.4) 


(A.5) 

V oo / 


-folds, A, is a monotonically increasing 


function of proper time t. However, its sign depends on the assumption of what the scale 
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ao represents. If oq is the initial value of a{t), such that a{t) > oq, then N starts from 
zero and take a large positive value when a{t) » ao- However, if one wants ao to represent 
the present value of the scalar factor, then N usually starts from a large negative number 
when a{t) <C ao and becomes positive only when a(t) > ao- 

Here we would like to present some exact solutions for some special cases, about which 
a general expression can be obtained for various parameters of the model, like, the field 
potential V(cr) and the coupling constant f(cr), in terms of the scalar field cr. 


8.1 x(N) = xo and h(N) = ho 

As a reasonable approximation, at late times, consider that the kinetic term K{a) oc 
so &/H ~ const. Specifically, when x = n^d) ~ ^Oi the field equations reduce to 

0 = —3 + Xo + y — ^{u — 2hu), 

0 = u" — {h + l)u + 2h{l + u — hz) — 2uh' + 2xo- (A.6) 

h = H/H^ is a kind of slow-roll variable, and thus may be treated as a constant in at least 
some region of field space. This is the case, for example, for power-law inflation, namely, 
a{t) oc where m < 1. For x = xq and h = ho, we find the following interesting solution: 

u{N) = uo + ui e^^°^ -|- U 2 (A.7) 

y{N) = yo + yie^^-^°^^, (A.8) 


where yo = 3 -|- ho{l — 5uo — houo), yi = 3^2(1 — 3ho), and uo,ui,U 2 are the integration 
constants. We also find that 


2x0 


-|- const, Xo = —ho{l + Uo — houo)- 


(A.9) 

(A.IO) 

(A.ll) 


Ka = TAi 

This further implies that 

(^0 + ill + 112 

Using (|A.9D , one may express these potentials as functions of the field a itself, namely 

U(ct) ~ + Vi (A.12) 

/(ct) ~ +/i + /2e“("/"°\ (A.13) 

where ^naoho = ^ combining the above expressions, we find 

A(cj) = ^ (a + B e^^°^ + C , (A.14) 
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where A = 3(l—uo)+ho(l—8uo—uoho), B = —3{l+ho)ui and C = —12hoU2- Equivalently, 


A(cj) = +Aie-2^/'"° (A.15) 

Thus, since ho < 0, it is the term proportional to or which is of greater 

interest at early {N < 0) or late {N » 0) times. In particular, when ho ~ 0, we find 

Via) = 3MlHlil + U 2 e^), VcBi^y) = ^mIhUuo + U 1 +U 2 e^). (A.16) 

In this case, since A(cj) ~ 3M|,i7g (1 — uo — ui), A(cj) acts as a cosmological constant term, 
which is positive for uq + ui < 1 . 


8.2 u{N) = Uo and xiN) = xq 

Let us recall that in the original action 

/? S 

^ ^ (2 + h), (A.17) 

‘iU‘2 

fi(y)Q = 2Afia)H'^ (l + h) = — 5 - (1 + h) u{a). (A.18) 

In a general situation, the coupling constant f{a) increases with logarithmic time, N, while 
the Hubble expansion rate, HiN), decreases with N. To this end, let us consider the case 
where the coupling constant f{a) scales as 1/H^, so that the function u(cj) is constant, i.e. 
u{a) ~ const = uo- In general, we would require uq < because any contribution coming 
from the higher order curvature corrections in low energy string effective actions should be, 
at least in a low energy scale, an order of magnitude smaller than the contribution from 
the Einstein-Hilbert term. Also, the function u(iT(iV)), which is likely to depend upon the 
gauge coupling strength, has to be an extremely slow varying function of proper time so 
that it assumes a nearly constant value. 

Let us further demand that the kinetic term iL(cj) scales as H^, so that x{N) = = 

xo- In this case, the explicit solution is given by 

yiN) = 3 — xo — QuohiN), hiN) = hi + ptanh. PiN — No), (A.19) 


where 


hi = 


Uo + 1 
2uo 


^ ^ V jup + 1 )^ + AxqUq 


2uo 


Further, a simple calculation shows that 


v /2 


Xo 


a = ±—. —N + const, 

K V 7 

H = cosh piN - Np), 


(A.20) 


(A.21) 

(A.22) 


where Np is arbitrary. For a canonical scalar, 7 > 0, we have xq > 0. A smaller value of 
Xo makes the potential flatter and hence increases the period of inflation. 
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Using ( [A.91 ), or equivalently N = 
to be Nq), we find 



(TK+ const (we will choose this last constant 


ua + 1 / 7 

^ 9,^9 — - J-rT^crn 

tZ A/tZ - 


V{a) = HlM^e “o V cosh""/3 ( \l 


—xo — 3uo — 6tto/3tanh/3 ( 4 / 


VV 2a:o 


/(o-) = 




e “0 V 2^0 sech^/3 


7 

2 x 0 ' 


-(TK 


(A.23) 

(A.24) 


Here a natural choice of uq is — 1 < uq < 0, so that the Hubble parameter H{a) is 
a monotonically decreasing function of the logarithmic time ln(a), or N. The effective 
potential is given by 


k^K{(t) = hImI e' 


“0 + 1 


'J ^^ 0 " 


' cosh^ f5 



X 


^3 — Xo — 3uo 


9uo 


hi + (3 tanh/3 



(A.25) 


In theories of quintessence, one Taylor expands the potential A((t) about the minimum, 
so as to obtain a dynamical scale for the mass of the quintessence field. This potential 
therefore merits further study. 


8.3 u{N) = uq and h{N) = /iq 

Consider, again, the case where u = uq, but now instead of x{N) ~ const, we demand that 
h{N) = H/H'^ ~ const. The solutions for x{N) and y{N) are given by 

X =-/io(l + uo - houo), y = 3 +/io(l - 5uo - houo). (A.26) 

One also finds the following relationships: 

a = ± — 4 [^N + const, H = . (A.27) 

K V 7 

The effective potential is then given by 

A(ct) = (3 + ho - 3uo - S/iqUo - hguo) ~ Aq e“'"/'"o. (A.28) 


8.4 x{N) = Xo and y{N) = yo 

Finally, as one more special case, let us consider that both the potential term V{a) and 
the kinetic term K{a) = scale with presumably, with different proportionality 
constants. For y = = yo and x = = xq, the solution is given by 


h{N) = 


3 + 5x0 - yo _ , 

— /T-q, 

Xo + 2/0 + 3 


M(a(lV)) 


(3 - Xq - yo)ho 
6 


— lii e 


‘IHqN 


(A.29) 
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where ui is an integration constant. The effective potential is therefore 


— T/„ o2/ioAf 


A(cj) = Voe 


3mi( 1 + ho)e2'*o^ + 2 





xo(a;o - 2yo - 3) \ 

5x0 - 2/0 + 3 J 


(A.30) 


From the above results, it is clear that at late times, a ^ cto, the leading order contribution 
to the potential comes from the term proportional to . 
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